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HAMILTONIAN OPERATORS AND ^*-COVERINGS 



P. KERSTEN, I. KRASIL'SHCHIK, AND A. VERBOVETSKY 

Abstract. An efficient method to construct Hamiltonian structures for non- 
linear evolution equations is described. It is based on the notions of variational 
Schouten bracket and £*-covering. The latter serves the role of the cotangent 
bundle in the category of nonlinear evolution PDEs. We first consider two 
illustrative examples (the KdV equation and the Boussinesq system) and re- 
construct for them the known Hamiltonian structures by our methods. For 
the coupled KdV-mKdV system, a new Hamiltonian structure is found and its 
uniqueness (in the class of polynomial (x, t)-independent structures) is proved. 
We also construct a nonlocal Hamiltonian structure for this system and prove 
its compatibility with the local one. 



Introduction 

We describe a method of constructing Hamiltonian structures for nonlinear evo- 
lution equations (or systems of such equations). The method is based on two 
concepts: the variational Schouten bracket and the £* -covering over a nonlinear 
PDE. 

In Section[l] we expose some general facts concerning the geometry of super PDE. 
In Section [21 we construct the variational Schouten bracket on a super version of 
Kupershmidt's cotangent bundle to a bundle and, following 2 , obtain an explicit 
formula for this bracket. In Section |31 simple computational formulas are deduced 
to check the Hamiltonianity of a bivector and compatibility of two Hamiltonian 
structures. Using the Schouten bracket, we define Hamiltonian evolution equations 
(including the cases when the Hamiltonian operator A may depend explicitly on 
time while the equation itself may not possess a Hamiltonian functional). This 
definition is equivalent to the operator equality 

where is the linearization of the equation and ^J. is the adjoint operator. To solve 
this equation, we introduce the notion of ^ ^.-covering (which is a particular case of 
a more general construction introduced in Sectional) and show that to any operator 
A satisfying the above equation there corresponds a function s on the £J,-covering 
such that £g{s) — 0, where £g is the lifting of the linearization operator for S' to 
the £J,-covering. In other words, the operators we are interested in are identified 
with shadows of nonlocal symmetries (of a special type) in the ^J,-covering. 

The reason to introduce the concept of £*-covering is two-fold. First, as it was 
just indicated, it allows to reduce construction of Hamiltonian structures to compu- 
tation of symmetries (with a subsequent check of additional conditions) for which 
a number of efhcient software packages exists. Second, to our opinion, this point of 
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view gives a new and fruitful insight into the theory of Hamiltonian structures for 
partial differential equations. 

In Sections [5] and [6l these methods are applied to the known examples of the 
KdV equation and the Boussinesq system. In Section [71 we construct a Hamil- 
tonian structure for the coupled KdV-mKdV system [B]. We also prove that this 
structure is unique in the class of (x, t)-independent polynomial structures. Never- 
theless, extending the initial setting with certain nonlocal variables, we find another 
Hamiltonian operator that serves a Hamiltonian structure for 'higher' coupled KdV- 
mKdV equations. This structure is compatible with the local one. It is to be noted 
that the theory of nonlocal Hamiltonian structures is not sufficiently developed yet 
and needs additional research. 

In the Appendix, we briefly recall the construction of the recursion operator 
for the coupled KdV-mKdV system (obtained earlier in [B]) by which the above 
mentioned Hamiltonian structures are related to each other. 

1. Generalities: Jet bundles and differential equations 

Let us formulate the main definitions and results we will use. For more details 
we refer to [2 [HI [2]. 

1.1. Jet bundles. Let tt: — > M be a vector bundle over an n-dimensional base 
manifold M and tToo : J°° (tt) M be the infinite jet bundle of local sections of the 
bundle tt. 

In coordinate language, if . . . , x„, u^, . . . , u™ are coordinates on E such that 
Xi are base coordinates and are fiber ones, then tToo : J°°{'^) — >■ M is an infinite- 
dimensional vector bundle with fiber coordinates m:^, where t = ii . . . i\-r\ is a sym- 
metric multi-index. 

Now, we generalize the definition of the jet bundle to the case of superbundles. 

Definition 1. Let E' be a supermanifold of superdimension {n + mo)\mi, and 
TT : E ^ M he a vector bundle over an n-dimensional even manifold M. If tt is split 
into the direct sum of two vector subbundles tt = tt" © tt^ such that the fibers of tt" 
are even and the fibers of ir^ are odd, then we say that tt (along with the splitting) 
is a superbundle. 

For a superbundle tt, we define the infinite jet superbundle tToo : J°°(7r) — > M by 
setting: 

(TToo)" - (7r")oo, 

where the superscript ^ denotes the reversion of parity. 

Denote by ,^{tt) the superalgebra of smooth functions on J°°{tt). 

Remark 1. By definition, we have 

-n^) = ^(tt") ®C3.(m) A*(^un((^^)")), 

where ^iin( • ) C J^( • ) is the subspace of functions linear along fibers. 
In what follows we shall use the term 'bundle' to mean 'superbundle'. 

1.2. Tiie Cartan distribution. Consider a bundle tt: E — ;> M and define the 
C°° {M)-supermodule r(7r) of its 'sections' as follows. If tt is even, then r(7r) is 
the module of sections of tt. If tt is a general superbundle, then we put r(7r) = 
r(^)0 ® r(7r)i, with r(7r)" = r(^0) and r(7r)i = (r((^i)n))n. 

Remark 2. Thus, in line with our definition of jets of superbundles, we define 
elements of r(7r) to be pairs of sections of 7r° and tt^. 
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Next, we note that every fiberwise linear function / on infinite jets J°°{'k) can 
be naturally identified with a linear diff'erential operator V/ : r(7r) C°°{M) and 
vice versa. Indeed, for even bundle tt the correspondence is given by the relation 



where s € r(7r), joo(s) is the infinite jet of s, a G M. The general case reduces to 
the even one, since =^iin(7r) = =^iin(7r° ® (tt^)^). 

Remark 3. The maps 



give rise to a map of supermodules joo '■ r(7r) — >■ r(7roo)- 

The infinite jet bundle tToo- J°°{tt) M admits a natural fiat connection such 
that the lift X of a vector field X on M is uniquely defined by the condition 



Vector fields of the form X generate an n-dimensional distribution on J°° (tt) called 
the Cartan distribution and denoted by ^(tt). Obviously, the Cartan distribution 

is Frobcnious in the sense that ["^(tt), "^^(tt)] C '^^(tt). In coordinate language, the 
Cartan distribution is spanned by the total derivatives. 

1.3. Horizontal calculus and evolutionary fields. Let ^: B — > M be a vector 
bundle and 7r^(0 : Bxm J°°(7r) -)■ J°°(7r) its puUback along tToo- The C°°(J°°(7r))- 

supcrmodulc r(7r^(^)) is defined as above: T{tt*^{£,)) = T{Cl (^c°-{m) C°°(J^(7r)), 
if ^ is even, and r(7r^(0) = r{n*^m° ®n^*oo{0)\ with r(7r^(0)° = r{n*^m 
and r(7r^(0)i = (r((7r^(0^)"))" if ^ is a general superbundle. 

Definition 2. A C°°(J°° (7r))-(super)module P of the form P = T{tt*^{^)) is said 
to be a horizontal module. 

Example 1 (horizontal forms). Let ^ be the qth exterior degree of the cotangent 
bundle to M. The corresponding horizontal module r(7r^(^)) is called the module 
of horizontal forms and is denoted by A'(7r). In coordinates, horizontal forms are 
generated by the forms / dxi^ A • • • A dxi^ , f G ^{i^)- 

Definition 3. Let P and Q be C°°(J°° (7r))-(super)modules. A map A: P Q is 
called ^ -differential operator (or horizontal operator) if it can be written as a sum 
of compositions of C°°(J°°(7r))-linear maps and vector fields of the form X. 

In coordinates, "^-differential operators are total derivatives operators. 

Example 2 (the horizontal de Rham complex). We define the first horizontal 

de Rham differential d: .^(tt) Ai(7r) = K^{M) ®c°°(M) by the formula 

d{f)(X) — X{f). In coordinates, we have d{f) = Di{f) dxi. 

The general horizontal differential d: A^(7r) — )■ A*+^(7r) is defined by the usual 
rules: 



V/(5)(a) = /(.7oo(s)(a)), 



ioo:r(7r0)^r(0, 



In coordinates, the lift of d/dxi is the ith total derivative 

d ^ d \ - , d 

— Ft. — L \ iiJ 




dod = 0, 

d{oJi A i02) = doJi Ai02 + (— l)^wi A doj2 
The differential J is a ^-differential operator. 



wi e A«(7r). 
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The cohomology of the horizontal de Rham complex 

^ ^(tt) 4 Ai(7r) ^ • • • a A"(7r) ^ 

are called horizontal cohomology and denoted by H'^{tt). From Vinogradov's 'one 
line theorem' [H [H [H] it follows that i79(7r) = i7«(Af) for g < n - 1. 

All "^-differential operators from P to Q form a C°°(J°° (7r))-(super)module de- 
noted by -rDiff (P,Q). 

Clearly, if P and Q are horizontal, then so is '^Diff(P, Q). 

Given a horizontal module P, let us define the horizontal infinite jet bundle 
Tip: J°°{P) J°°(7r) as follows. If P is even, then the fiber of J°°(P) over 
9 e J°°{-k) consists of equivalence classes, denoted by j{p){0)^ of elements p £ P. 
Two elements pi and p2 are equivalent if their total derivatives of all orders coincide 
at 9. For a general horizontal supermodule P, we as always define iTp — npo and 
TTp = TT^i^n- Correspondingly, r(7rp) — r(7rpo) r(7r(pi)n)'^. 

Clearly, the horizontal jet bundle vrp: J°°{P) J°°{n), P = r(7r^(0), is 
isomorphic to the pullback 7r^(Ccx)): J'^iO J°°{tt) — > J°°{tt) and, thus, r(7rp) 
is a horizontal module. 

Similarly to Remark [31 we have the natural operator joo - P J°° (P) ■ 

For every ^-differential operator A : P — Q there exists a unique homomorphism 
of C°°(J°°(7r))-supermodules Ha: J°°(P) J°°iQ) such that the diagram 



P 



is commutative. 

Let us recall the definition of adjoint operator. Consider A e "^Diff (Pi, P2). 
The adjoint operator A* G "TDiff (P2, A), P = Hom^(^)(P, A"(7r)), is uniquely 
defined by the equalitjQ 

(p,A(p)) = (-i)^^^(A*(p),p), pgA, pePu 

where (• , •) is the natural pairing P x P ^ H"{tt). 
In coordinates, we have 



(1) 



where a^^ G '^(t'')? the superscript 'st' denotes the supertransposition, and Dr ~ 
Ai o ■ ■ • o for r = ii . . 

Equivalently, adjoint operator can be defined using the following fact. Consider 
a horizontal module P and the natural complex 

'rDiff(P,^(7r)) -> 'g'Diff(P, Ai(7r)) ^ "^Diff (P, A2(7r)) ^ ■ • • 

^ '^Diff(P, A"-i(7r)) ^ <rDiff(P, A"(7r)) ^ 

with the differential A 1—^ do A. Denote its cohomology by iJ^(P). We have 

(hi{P)^0 forO<g<n 

|i7"(P) = P. ^ ^ 

Each '^^'-differential operator A: P — Q gives rise to a cochain map between two 
such complexes. The corresponding map of the nth cohomology A* : Q ~^ ^ is the 



"'^Here and below, symbols used at the exponents of (—1) stand for the corresponding parity. 
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adjoint operator. Note that the natural projection fi: "^^Diff (P, A"(7r)) P has 
the form A ^ A*(l). 

Recall the most important properties of adjoint operators: 

(1) A and A* are of equal parity; 

(2) (A10A2)* = (-l)^i^^A^o^At; 

(3) A** = A (here we identify P and P). 

A vector field Z on J°°{7t) is called vertical if ^1,^00(4,/) = 0. For a horizon- 
tal module P a vertical field Z generates a natural action Z: P ^ P, which in 
coordinates is the component-wise action. 

A vertical vector field Z is said to be evolutionary if [Z, X] = for all vector 
fields X on M. 

It is easy to see that evolutionary fields are uniquely determined by their re- 
strictions to ^\in{E), where E is the space of the bundle tt: E ^ M. Moreover, 
the map Z i— > Z\^^_ is a bijection between the set of all evolutionary fields 
and Homc=c(M)(=^iin(i^),=^{7r)). We identify Y{.om.c^{M){'^\\n{E), ^{tt)) with the 
horizontal module r(7r^(7r)) and denote it by k{'k). 

In coordinate language, the evolutionary field that corresponds to a vector func- 
tion if — {ip^ , . . . , ip™) has the form 

J.T 

Let P be a horizontal module. The linearization of an element P e P is a 
"^-differential operator Ip: >t{Ti) P defined by the formula 

Denote by the square brackets the horizontal cohomology class of a horizontal 
form. Since evolutionary fields commute with the horizontal differential, the coho- 
mology class [9^{ijj)\ for uj e A"(7r) is well defined by [w]; denote it by d^{[uj]). 
By ((TJ we have 

where £: A"(7r) >c{t^), £(w) = is the Euler operator, which takes La- 

grangians to the corresponding Euler-Lagrange equations. Of course, the value 
£(a;) is completely determined by the cohomology class [w]. 

In coordinates, E{Ldx^ A • • • A da;") — {SL/Su^, . . . ,SL/du"^), where 5L/5u^ = 

E.(-l)l^l^r(9V5<). 

Remark 4. From Vinogradov's 'one-line theorem' [12l[l3l[T4] it follows that 

(1) ker£/J(A"-i(7r)) i7"(M); 

(2) G im£ if and only if £^ ~ t^p- 

1.4. DifTerential equations. Again, consider an element P of a horizontal mod- 
ule P. The locus 

'^°° = {joo(F)=0}c J°°(7r) 

is called differential equation defined by P. We assume that the natural map 
^ M is a subbundle of the bundle tt: J°°(7r) M. The restriction of the 
Cartan distribution to is denoted by "^((f ). Clearly, dini^(<o') = dim "if = n. 

Example 3 (evolution equations). Consider the bundle tt : P x M — s- A/ x M. Denote 
the coordinate along M by i. Then 9$ — Df — d/dt is a canonical evolutionary field 
on J°°(7r). In coordinates, $ = {ul, . . . ,m™). Let 9<p(t) be a family of evolutionary 
fields on J°°(7r). The equation C J°°(7f) given by the element P = $ - ip{t) G 
>f(7f) is called evolution equation. In coordinates, it has the form ut = </'(t). Note 
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that = J°°{'k) X M, with the Cartan distribution generated by that on J°°(7r) 
and the field Dt = dldt+ 9ip{t)- 

The restriction of the hnearization ^i? to the equation S°° is called the lineariza- 
tion of and is denoted hy £s- x P, where x = x{Tr)\^^. 

An evolutionary field tangent to (S""" is said to be a symmetry of the equation. 
Obviously, d^p is a symmetry if and only if igi^) = 0, G >f. 

The horizontal de Rham complex on J°°(7r) can be restricted to (0°°. Its coho- 
mology are called horizontal cohomology of equation S"^ and denoted by H'^ (£') . El- 
ements oi H"^^^ {S") / H"^^ (M) are conservation laws oiS"^. If the equation at hand 
satisfies the conditions of Vinogradov's 'two-line theorem' [12l [131 E] then there is 
an inclusion i7"-i(^)/i7"-i(M) ^ ker£J. The element i{r]) e kertg C P that 
corresponds to conservation law 77 is called its generating function. 

In particular, evolution equations satisfy the conditions of the two-line theorem. 
In this case, i(r?) = £(770), where = ryo + Vi ^ dt, 770 G A""i(^), 771 G A""2(^). 
Thus, to find conservation laws of an evolution equation one has to solve the equa- 
tion tg{'4') = and choose those solutions -0 that fulfill the condition = 

Let S"^ and be two differential equations. A surjective map t: § ^ is 
called covering if it preserves the Cartan distribution. 

Example 4. A horizontal jet bundle ttp : J°°{P) — > J°°(7r) is a covering. A gen- 
eralization of this example will be discussed in Section |4l 

Example 5. Let (0°° be given by an element F. Consider equation S""" 

\dr/dxi=gi, 
[dr/dxn = gn, 

where (71, . . . , are functions of xi, . . . , x„ and u^^. If the compatibility condition 
A(<?j) = D,{g,) on ^ d7^\^^ =0 7? = 5. dx, £ A^,^) (4) 

holds true, then the natural projection r: S'°° — S'°° is epimorphic. Obviously, 
T preserves the Cartan distribution, so that it is a covering. Thus, each closed 
horizontal 1-form gives rise to a covering of the form ([3]). In particular, when 
n = 2, condition ^ means that 77 represents a conservation law of (0°° . The new 
dependent variable r is called nonlocal variable. 

In a similar way, we can define a covering over (0°° corresponding to a closed 
one-form on this equation, etc. In this manner we construct particular coverings in 
Sections EMU 

Clearly, each "^-differential operator A on (0"°° can be lifted to a "^-differential 
operator A on (0°° . In particular, we have the operator is on A symmetry 

of <?°° is called a nonlocal symmetry of S""" in the covering under consideration. 
Solutions of the equation (si^^) — are called shadows of nonlocal symmetries 
of (0°° in this covering. In a similar way, since the horizontal de Rham differential 
is a ^-differential operator, we can lift the horizontal de Rham complex to S'°° and 
construct the theory of nonlocal conservation laws in our covering. Solutions of the 
equation i*g{i}}) — are called nonlocal generating functions. 

2. Variational Schouten bracket 

We start with a super version of Kupershmidt's cotangent bundle to a vector 
bundle [S]. 
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For a vector bundle tt: E — > Af, dimAf = n, we consider the bundle tt: E ~ 
E* (g)M A"(T*M) ^ M, where £;* ^ M is the dual bundle to E M, and the 
superbundle J^: J^'' = tt (even subbundle), J^'^ = tt (odd subbundle). 

The superbundle Xo : J°°(jr) ^ M 




M 

is called the cotangent bundle of the bundle tt. It is clear that 

J-(^) = J-(r«(7roo))). 

Denote by , j = 1, . . . ,m, the fiber coordinates in E dual to with respect 
to a volume form on M (they are sometimes called 'antifields'). Then Xi, m^, 
will be the coordinates in J°°(J^), with z^, being even and being odd. 

It is clear that >!:(J^)° = xje- and x{,^Y — where >ix = T{J(^*{n)). 

Define the variational Schouten bracket (antibracket) on the space H"{J^) by 
putting 

lF,Hl = {E{H),a{E{F))), F, H e H^{.je), (5) 

where £ is the Euler operator and the operator a : k{J(f) — >c{J(f) acts according 
to the formula a(V', = (Vj for ip G xje and -0 £ >(j(r- In coordinates, we get 

It is readily seen that the variational Schouten bracket defines a Lie superalgebra 
structure on H"{J(f): 

[F,Fl = -(-l)(^+i)(«+i)lF,Fl, 

Remark 5. A different concept of the Schouten bracket (acting on a different space) 
the reader can find in 15, p. 226]. 

Denote by '^DiH^^^^ {P, Q) the module of fc-linear skew-symmetric '^-difi'erential 

operatorsFx- • -xP Q. The subset '^Diff(t)"'''^(P, P) C "^Diff (t°'^(P, P) consists 
of skew-adjoint in each argument operators. 
Let us define multiplication 

'rDifff,^)''"(P,^(7r)) X <^Difff;^<=-(P,^(7r)) ^ '^Diff(t;7)(P,^(^)) 

by setting 

(AiA2)(pi,...,_Pfc+;) = X! (~l)'^^l(P<T(l,fc))A2(p„(fe+i,fc+i)), 

where Ai G 'g'Diff (t)™(P, ^(tt)), A2 G -^Diff (",'^''^(P, ^(tt)), 5^ C Sn is the set of 
all (i, n — i)-unshufHes ([in])j i-e., all permutations a G Sn such that a{l) < a{2) < 
■ ■ ■ < a{i) and <7{i + 1) < cr(i -|- 2) < • • • < cr(n), (—1)'^ is the sign of permutation cr, 
and Pa{kiM) stands for p„(^ki),- ■ ■,Pa(k2)- 
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Next, since by definition elements of --^(tt) are identified with differential oper- 
ators from r(7r) to C°°{M), we have the natural inclusion "^Diff (i<-(7r), ^(vr)) 
^(J^), which uniquely prolongs to the isomorphism of algebras 

'^Diff(tp(i<(7r),^(7r)) ^ ^(Jf). 

Using ([2]) we can show in a standard way that 

= ^^Difffi')^'^(i<(^), ^(n)) © (6) 

Below, we use the shorthand notation >c — >c{tt). 

Now, following [5] , we want to compute the variational Schouten bracket in terms 
of skew-adjoint "^-differential operators. 

To this end, note that from the definition of the Euler operator it follows that 
its restriction 

e|^Diff-w(^.,A"W) : <^Difff,^p(iV,A"W)^'^^Diff(';^p(i.,>.)©'^^Difff^T)(^,>^), 

has the form £|^Diff =,''™(ii,A"(7r)) (^) = ('7(A), (-l)*^" V(A)), where 

77(A)(Vi,...,^fc)=^l.^,,...,^,(l), 
^i{A){^u . . . , V-fc-i) = (A(^i, . . . , ^fc-i))*(l), 

ipi e >(, iA.i,i,...,i,k{v) = 9^{A){iJi,...,ipk)- 

In coordinates, 77 = ((5/'5m\...,(5/'5w'") and = {-l)''-\S/5p\ ... ,5/6p'^). 
We can rewrite 77 in the following form: 

77(A) = 77(m(A)), a e <^Diff(tp(>., A"W), 

where 

77(n)(7^i, . . . ,7Afe) = ih^^,_^,_,{^Pk), □ G ^Bmit%iK,M). 

Indeed, take the equality 

[A(7Ai,...,Vfc)] = (□(^i,...,^fc-i),7/'fe), D^MA), 
and apply 9^ to both sides. This yields 

[9^(A)(7/>i,...,7/>fe)] = (9^(n)(V^i,...,7/;fe_i),7Afe), 

and so 

(¥',^k^„...,^,(l)) = (^,^a^,....,^,_,(V'fe))- 

Thus, for the variational Schouten bracket of two operators A G "^Diff ^'^^i^^ {>c, k) 
and B e '^#'Diff(f_:i^(>f, x) we have 

(|A,Bl(l/.l,...,7Afc+i-2),V'fc+i-l) 

= (-1)'"' E (-l)'(^B.VMi.-i,(^-(0)'^(^-('+i,fe+'-i))) 

^ (-l)'^(£lvva,.-i,(^-w)'^(^-('=+i^'=+'-i)))- 
Here and below we assume that 5*,^ = if i < or i > 77. 

Let us split the sums obtained into two parts depending on whether a{k + l — 1) = 
fc + / — 1 or not: 

(|A,5l(7/.l,...,7Afc+i-2),V'fc+i-l) 
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Thus, we have 

[ABi(^i,...,Vfe+i-2)- J2 (-ir^B,^.,,,,_„(A(^.(i,fc+i-2))) 

+ E (-i)"A(^B.v,„<,,,_,,(^.(o)''^^(;+i,fc+;~2))- (7) 

From the definition it immediately follows that 

{A, Wl (V'l , . . . , ?Afe_2) = A(£ (w) , , . . . , 

for uj G H^{tt); in particular, |iy9, w] = ((^, £(w)) = 

3. HAMILTONIAN EVOLUTION EQUATIONS 

An operator A G "^Diff (iV, x) is called Hamiltonian if |A, A] = 0. As in the clas- 
sical Hamiltonian formalism, a Hamiltonian operator defines a Lie algebra structure 
on //"(tt) via the Poisson bracket 

{c^i,c^2}a-(A(£K)),£(cj2)). 

Remark 6. Hamiltonian operators are uniquely determined by the corresponding 
Poisson brackets. 

Remark 7. A Hamiltonian operator A gives rise to a complex 

^ iJ"(7r) ^ ^ -^^Diff ^ '^Diff("2^)-'^^(iV, >f) ^ ■ • • , (8) 

where (A) = |A, A], called the Hamiltonian complex. 

Formula ([7]) yields a well-known criterion for checking a skew-adjoint operator 
to be Hamiltonian (see, e.g., [1113): 

[A, AKV^l, V2) = -^A>i(A(^2)) +^A,^2(^(V'l)) -^(^1^,(^-2)) = 0. 

Another practical way to check the Hamiltonian property of an operator is to 
use formula ([5]). In coordinates, it gives: 

E SWa SWa _ Q modulo total derivatives (9) 

3 

„ , ^ .5Wa SWa 
or t 



Sui dpi 
3 
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where Wa G {,J^) is the element that corresponds to the operator 2 A under the 
isomorphism ([5]). In coordinates, the element Wa for an operator J^t'^t^t has 
the form Wa = Z^i.j.r Pip'- 

The condition for two Hamiltonian operators A and i? to be a Hamiltonian pair, 
i.e., [A 51 =0, is 

^^^swASWB^sn^swA\ 



^-^ V Sui 5pi 5ui Sp- 
Note that the skew-adjointness in terms of Wa amounts to the equality 

E^p--™'- (10) 

3 

Let yl be a Hamiltonian operator. Evolution equation ut = f is said to be 
Hamiltonian with respect to A if 

At-IA/l-O, (11) 

where At = dA/dt (both A and / can depend on the parameter t). 

If A does not depend on i, then for each H E _ff"(7r) the evolution equation 
Ut — A{£,{H)) is a Hamiltonian evolution equation. The element H £ i/"(7r) is 
called the Hamiltonian. Notice that in this case condition (ITU) means that / is a 
1-cocycle in the Hamiltonian complex (|5]). A Hamiltonian H exists if and only if / 
is a coboundary. 

If a Hamiltonian H exists and does not depend on t, then we have 

Dt{H) = 9fiH) = 9a(e(h)){H) = {H,H}a = 0. 

Thus, there exists a conservation law given by 77o+?7iA(ii, 770 G A"(7r), rji G A"^^(7r), 
such that [770] = H E H"{tt), where [rjo] is the cohomology class of 770 in _ff"(7r). 
In other words, the generating function of this conservation law equals £-{H). This 
conservation law is called the conservation law of energy. 

Theorem 1. Let he an evolution equation Ut — f which is Hamiltonian with 
respect to a Hamiltonian operator A. Then we have 

£soA + Aotg = {). (12) 

Proof. By ©, 

{At - [A, fim = Ati^p) + iAM - Mtfm - ^Mm, 

thus 

At-{AJ\=At + 9f{A)-Aotf-tfoA. 
Hence, {d/dt + 9f - tf)o A~ Ao {d/dt + 9/ + £}) = 0. It remains to note that 
£s = Dt-£f=d/dt+9f-£f. □ 



Remark 8. For equations possessing a Hamiltonian, relation (|12p can be found 
elsewhere (see, e.g. [H]). 

We call solutions of (I12p variational bivectors on the equation under consider- 
ation; Hamiltonian operators that make a given equation Hamiltonian are, thus, 
special variational bivectors on the equation. Obviously, variational bivectors (and, 
in particular, Hamiltonian operators) take generating functions of conservation laws 
of the equation at hand to symmetries of this equation. 

Proposition 1. Let S'°° be an evolution equation Ut — f . If two operators A, 
A' e "^Diff (i<:, x) satisfy the equation 

i^oA + A'otg^O, (13) 

then A! = A. 
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Proof. Rewrite ([T^ in the form 

{Dt-ef)oA-A'o{Dt + e})^o. 

Commute the right-hand side of this equahty with the operator of multiphcation 
by t. This gives A' = v4. □ 

4. A-COVERINGS 

In this section, we describe a construction, that reduces solution of equation (jl2p 
to finding shadows of nonlocal symmetries in a special covering over (o (the ^ J- 
covering). 

Let S'°° be a differential equation, and A: P Q he a. '^^ -differential operator 
between two horizontal modules P and Q over (0°°. Consider the homomorphism 
/ia: J°°(P) J°°(Q) that corresponds to A. If Ka = ker/iA C J°°{P) is a 
subbundle of J°°{P), then = 7rp|^,^ : Ka ^ is a covering. We call it 
IS.- covering. 

In terms of local coordinates, if A = || X)r '^t ^t|| and are fiber coordinates 
of a, where a is such that P — T{a)^ then A-covering is defined by the equations 

Y,(^rWi=^- (14) 

We can think of fibers of A-covering as even or odd. Here we prefer the latter 
viewpoint, so that ^a is a superbundle. 

A-coverings are useful mainly due to the following obvious fact. 

Proposition 2. Let R — r(7) be a horizontal module over S"^ . Then there is an 
isomorphism 

riin(fcl(7)) = "^BiUiP, R) / {V £ <^Diff (P,i?) I F = Do A,n e ■^Diff(g,i?)}, 
where denotes space of fiberwise linear sections. 

Proof. The isomorphism takes s G r\in{k'^{'j)) to the equivalence class of the oper- 
ator Vs'. P ^ R given by the formula 

Vs{p) = SOJ^{p), 

where s is an extension of s to J°°{P), p <E P. □ 

In coordinate language, Dr at the jth component of the operator goes to w^. 

Now suppose that we are given a "^-differential operator V : R ^ R' over S"^ . 
Let us lift it on Ka and consider ker V. In view of Proposition [21 we can identify 
fiberwise linear elements of ker V with solutions V G "^Diff (P, P)/{ □ o A } of the 
equation 

VoV^V'oA (15) 

Thus, (jlSp amounts to the equation V(s) = on the A-covering. 

In particular, equation (IT3l) is equivalent to the equation dsi'f) = on the £J- 
covering, where if is fiberwise linear vector function. Note that in this case the 
^J-covering can be identified with the cotangent bundle J°°{J(f) x M. Under this 
identification, the Cartan planes on K^^ are spans of the Cartan planes on J°°{,J^) 

and Dt = d/dt -\- dp where / — {f,£j{'w)), if the equation at hand is Ut — f. 
Moreover, the equivalence classes of operators from Proposition [5] are in one-to-one 
correspondence with "^-differential operators on J°°{tt). 

Remark 9. From the above said, we see that Hamiltonian operators are shadows of 
nonlocal symmetries in the ^ J-covering. 
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Remark 10. Solutions V E "^^Diff (P, iT!)/{ □ o A } of equation ^ can be found 
straightforwardly and the computations will be essentially the same as when one 
solves the equation V(s) =0 on the A-covering. Nevertheless, in our computations 
we prefer the second approach, because in the case V = it reduces the problem 
to finding shadows of nonlocal symmetries (see above) for which efficient software 
exists. 



5. The Korteweg-de Vries equation 

Here we show how the above introduced techniques work with a simple and 
well-known example of the KdV equation 

Ut = Uxxx+UUx- (16) 

Local coordinates in S'°° are 

X, t, U = Uq, . . . ,Uk, ■ . ■ , 

where Uk = d^u/dx^ (similar notation is used in the subsequent sections as well). 
In these coordinates, the total derivatives are 

fc>0 
fc>0 

The £g-covering. The linearization operator for (|16l) is 

(s^Dt-Dl-uDx-ui, 
while the adjoint is expressed by the formula 

-Dt + Dl + uDx. 

Following the general scheme, we construct the £J,-covering by introducing the 
odd variables p = po, pk — D^{p) that satisfy the equation 

Pt =P3 + upi. 

Solving the defining equation. Let us now extend the total derivatives up to the 
total derivatives on the £J-covering 

d 

Dx = Dx + ^Pk+i-^, 

Dt^Dt + y^Dlip^+upi)^. 

Then, solving the equation lg{F) = 0, that is, 

bt{F)=bl{F)+ubx{F)+uiF, (17) 

with respect to the function F — -Fipi, where Fi — Fi(x,t,u, . . . ,Uk), we obtain 
two independent solutions 



2 1 

F°=pi, F^ = p3 + -upi + -uipo 



to which there correspond two "^-differential operators 

A^'^Dx, ^Dl+'^uDx + ^'^ 

the classical Hamiltonian structures for the KdV equation 
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The Hamiltonianity test. To demonstrate how the method works, we shaU check 
the Hamiltonianity of the operators and in a straightforward way. Obviously, 
both operators are skew-adjoint. 

For A^ , the corresponding bivector Wq = Wa" is 

Since SWo/du = 0, we get 



\ Su Sp / 



= 0. 



For A^, one has 



Wi = FVo = {P3 + gWPi + gUiPojpo = PsPo + gWPlPO- 



Consequently, 



and 



6Wi 2 
du 6 



SWi _dWi ^ dWi ^^dWi ( ,2 \ ^ (2 
5p dpo dpi ^ dpz 



(p3 + 3"Pi) - D^(^-upQ^ - DKpq) 



4 2 
-2p3 - ^upi - gWiPo- 



Hence, 

JVFi SWi _4: ^ /4 

(5u (5p 3 
that implies £.{SWi/5u ■ 6Wi/Sp) = 



PqPiP3 = Dx(^-^PqPiP2^ 



Remark 11. In ^ we describe a class of equations which have the property that 
([T^ automatically implies the Hamiltonianity. In particular, KdV belongs to that 
class, thus the above verification might be skipped. 

Nonlocal Hamiltonian structure. Let us introduce a new (odd) nonlocal variable 
determined by the equations 

rx = uipo, 

rt = UiP2 - U2P1 + {uui + U3)po 
(see Example O. Then an additional solution of equation PT|) arises: 

2 4 /4 2 4 \ /4 1 \ 1 

P =P5 + -^up3 + 2miP2 + l^gU + -U2 jpi + [gUUi + -U3 jpo - -uir, 

to which there corresponds the operator 

A^ = Dl + ^uDl + 2mDl + (^^2 ^ iu2)D^ + i^^um + i^a) - \uiD-^ o m. 

Remark 12. Here and below we use the following correspondence between nonlocal 
variables and operators (in the case of evolution equations with one-dimensional x) . 
Let be the variables in the fibers of the £J-covering and nonlocal variable r be 
determined by the relations 

Ux = Ek^jalpl, 

(cf. Example [SJ. Then the corresponding operator acts on ip — {ip^, . . . , ip"^) by 
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Is 



Simulating the techniques developed for the local case, it is a straightforward 
check that is a Hamiltonian structure and all three structures are pair-wise 
compatible. Moreover, they are related to each other by the classical recursion 
operator 

i.e., — Ro AS' and A? — Ro . In a similar way, one can find a whole infinite 
series of nonlocal Hamiltonian structures for the Korteweg-de Vries equation. 

Remark 13. We stress here the word simulating above: at this moment, we do not 
have a consistent theory of Hamiltonian structures in the nonlocal setting. We hope 
to develop it elsewhere. 

6. The Boussinesq equation 

In this section, we shall present, as another illustration of the above developed 
methods, computation of local and nonlocal Hamiltonian structures for the classical 
Boussinesq equation. We consider this equation as the system of the form 

Ut = UxV + UVx + CrVxxx, 

Vt=Ux+ VVx, (18) 

where cr G M is a constant. 

All computations presented below were done by the software system described 
in [71 Ch. VIII] and we expose here final results only. 

The fg-covering. The linearization operator restricted to is 

'^vDx - Dt crDl + Ml 
Dx vDx +vi - Dt^ 

while the adjoint one is expressed by 

-vDx ~ vi+ Dt ~Dx 
-aDl + Ml -vDx + Dt^ 

Hence, the ^J,-covering with the odd nonlocal variables pi, qi is defined by 

Pt = vpi + vip + qi, 

qt = o-p3 ~uip + vqi. 

Local Hamiltonian operators. In a completely similar way as described for the KdV 
equation in the previous section, we solved the symmetry equation in the ig- 
covering of the classical Boussinesq equation. We found three local solutions of 
the form 

= qu 
G'^Pi; 

= 2o-p3 + 2upi + uipo + vqi, 
= vpi + vipo + 2qi; 

— Aavp3 + 6avip2 + 2{3av2 + 2uv)pi + 2{av3 + uvi + uiv)po 
+ 4(753 + (4u + M^)gi + 2uiqo, 

= 4crp3 + (4u + v'^)pi + 2(ui + vvi)po + Avqi + 2viqo. 

In classical operator notation, they are represented as 

Ai^f^ .2 ^ {^<^Dl + 2uDx + ui vD, 

\Dx J' \ vDx + vi 2D, 
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while the third operator has the entries 

= 4:avDl + daviDl + 2{3aV2 + 2uv)D^ + 2(cru3 + uvi + Uiv), 
= 4:aDl + (4w + v^)D^ + 2ui, 
Al^ = 4aDl + (4u + v'^)D^ + 2(ui + vvi), 
= "ivDx + 2vi. 

Hamiltonianity and compatibility. To test Hamiltonianity and compatibility condi- 
tions for the operators A^, A^, A^, we construct the bivectors 

Wa^ = qiPo+Piqo, 

Wa^ = (2tTp3 + 2upY + MiPo + w9i)Po + (vpi + uipo + 2gi)go, 
Wa^ = (4crfp3 + 6(TWiP2 + 2(3(Ti;2 + 2uv)p\ + 2(av-i + uv\ + u\v)pq 
+ 4(793 + (4m + + 2?iiqo)po 

+ (4crp3 + (4u + v^)px + 2(ui + i't;i)po + 4i'qi + 2uiqo)go 
and straightforwardly check that 

{Wa^.WaA^^. l<i<j<3, 

i.e., the operators A^, meet both Hamiltonianity and compatibility condi- 

tions. 

Nonlocal Hamiltonian operators. In order to describe nonlocal results we introduce 
three new nonlocal variables ri, r2, rs over the ^J-covering by the following defini- 
tions 

ri,x =PoUi +qoVi, 

r\,t = P2<yvi - piav2 +Po{o-vs + uvi + uiv) + qo{ui + vvi); 
r2,x = Po(ctw3 -I- uvi +uiv) + qo{ui +vvi), 
r2,t = P2<y{ui + vvi) - pia{u2 + VV2 + vl) 

+ Po{(^U3 + 2avv3 + SaviV2 + uui + 2uvv\ + uiv^) 

+ 9o(o-W3 + uvi + 2uiv + v'^vi); 

f3,x = Po(4o'U3 -|- Qavvz + 12crwit;2 + 6uui + 6uvvi + Suiv'^) 
+ qo{^(TV3 + Guvi + 6uiv + Sv'^vi), 

'f3,t = P2CT(4crW3 -|- duvi + 6uiV + 3v^Vi) 

+ Pia{—4av4, — 6uv2 — 12uiWi — 6U2V — 3v^V2 — 6vvi) 

+ Po{^cr^V5 + lOauvs + 18auiV2 + 18o-U2Wi -|- lOau^v + 9av^V3 

+ 30(TTOiW2 + Q(Tvl + 6u'^vi + 12uuiv + 9uv'^vi + 3uiv^) 

+ qo{4:aU3 + lOavvs + 12cruiW2 + Quui + 12uvvi + 9uiv^ + 3v^vi). 

Using these nonlocal variables, we derived the following three nonlocal Hamiltonian 
structures given by 

= 8a^p5 + 2a{8u + 3u^)p3 -|- 6<t(4ui -|- 3vvi)p2 

+ 2{8au2 + 9avv2 + Qav^ + 4u^ + 3uv'^)pi 

+ (40-^3 -|- 6avv3 + 12av\V2 + 8uui + &uvv\ + 3uiV^)po + 12avq3 
+ 20aviq2 + {16aV2 + 12uv + v^)qi + 2{2aV3 + 2uvi + 3uiv)qo 
- 2uiri, 
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= 12avp3 + 16aviP2 + (12ctw2 + 12uv + v^)pi 
+ {Aavs + 8uvi + 6uiv + 3v'^vi)po + Saq^ 
+ 2(4m + 3v^)qi + 2(2ui + 3vvi)qo - 2viri; 

= 32cr2?;p5 + SOa^vipi + 8(t(14otj2 + Suv + v^)p3 

+ 4ct(22ctw3 + 24mwi + 24miv + 9v'^vi)p2 + 4(10cr^W4 + 20(TW2 

+ 26cruiVi + 16aU2V + 9av^V2 + 12avvl + Su^v + 2uv^)pi 

+ 4{2cr^V5 + 6cruv3 + llauiV2 + 9cru2Vi + iauav + Sav^va 

+ 12crvviV2 + 3avl + iu^vi + Suuiv + Suv'^vi + uiv^)po + IGa^qs 

+ 8c7(4u + 3v'^)q3 + 16cr(3ui + 5vvi)q2 + {32aU2 + 64avv2 + Uav^ 

+ 16u^ + 24uv'^ + v^)qi + 4(2crU3 + iavvs + 6aviV2 + 4uui 

+ 4uvvi + 3uiv^)qo — 4uir2 — 4(cru3 + uvi + uiv)ri, 

= 16(T^p5 + 8(7(4u + 3v'^)p3 + 16c7(3ui + 4vvi)p2 

+ {32aU2 + 4Savv2 + 28av^ + IGm^ + 2Auv'^ + v^)pi 

+ 4(2o-U3 + 4:avv3 + 8aviV2 + iuui + 8uvvi + 3uiv'^ + v^Vi)po 

+ 32(Tvq3 + A8aviq2 + 8(4crw2 + 4mu + v^)qi 

+ 4(2crw3 + Auvi + 4:Uiv + 3v^vi)qo — 4:Vir2 — 4(mi + vvi)ri; 

= -32ctV - 16cr2(6u + 5w2)p5 - 80cr2(3ui + 5oti)p4 

- 2cr(160cr?i2 + 280(jvv2 + 2M(Jvl + 48?!^ + 8Quv'^ + 5v'^)p3 

- 4o-(60o-U3 + llQavvs + 216crt;iW2 + 72uui + 120uoti + 60uiv'^ 
+ 15v^vi)p2 - 2(48<7'^U4 + lOOcr VV4 + 244(7 V1V3 + 168(7 V2 

+ 96(7UU2 + 200auvv2 + 136auvf + 68aui + 260auivvi + 80au2V^ 
+ 30av^V2 + GOav'^vl + 16u^ + AOu^v"^ + 5uv^)pi 

- (16c7^U5 + 40a'^vv5 + 120(7^vif;4 + 208cr^t;2U3 + 48auu3 

+ 120cruvv3 + 232auviV2 + 88cruiU2 + 220cruivv2 + 156c7Uii)f 

+ 180c7M2t^vi + 40c7M3t;^ + 20c7t;^t;3 + 120(7t;^t;it;2 + QOavvl 

+ 48u^«i + 80u^vvi + 80uuiv'^ + 20uv^vi + 5uiv^)po - SOa'^vq^ 

- 224:a^viqi - 40c7(8c7t;2 + 4uv + v^)q3 - 8a{30aV3 + 32uvi 
+ 30uiv + 25v'^vi)q2 - (96c7^t;4 + 192c7Mt;2 + 256(7Uii)i 

+ 160(7U2W + 160av'^V2 + 220avvi + 80u^v + 40uv^ + v^)qi 

- 4(4(7^z;5 + 12auv3 + 22auiV2 + 18(7M2Wi + lOausv + 10av^V3 
+ 30(7wwiU2 + 6<7vl + 8u'^vi + 20uuiV + lOuv^Vi + 5u-iV^)qo 

+ 2uir3 + 8{av3 + uvi + Uiv)r2 

+ 2(4(7M3 + 6(77;w3 + 12(7Wi?;2 + &uui + Quvvi + 3?iiw^)ri, 
= ~80a^vp5 - llQa^vipi - 8(7(28(7i;2 + 2{)uv - hv^)p3 

- 16(7(ll(7t;3 + lAuvx + lbui_v + lQv'^v\)p2 

- {80a^V4, + 160(7UW2 + 224(T«it;i + l&QaU2V + 12Qav^V2 
+ lAQavvl + 8Qu^v + AQuv^ + v^)pi - (16c7^t;5 + 48crut;3 
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+ 88auiV2 + 88fTU2Wi + 40fTU3W + AOav^vs + 160avviV2 + SGavf 
+ ASu^vi + SOuuiv + SOuv^vi + 20uiv^ + 5v'^vi)po - 32o-^g5 

- 16(t(4u + 5v'^)q3 - 48a{2ui + 5vvi)q2 

- 2(32c7M2 + S0(TVV2 + 52avl + 16u^ + 40uv'^ + 5v^)qi 

- 4(4(TU3 + lOcrvva + 16crviV2 + 8uui + 20uvvi + lOuiv^ 
+ 5v^vi)qo + 2vir3 + 8(ui + vvi)r2 

+ 2{4crv3 + 6uvi + 6uiv + 3v'^vi)ri. 

In the conventional notation the operators A^, A^, and A^ have the following 
entries: 

Af^ = 8a'^Dl + 2o-(8u + 3v'^)Dl + 6(7(4ui + 3vvi)Dl 
+ 2{8aU2 + 9avv2 + 6avf + Au^ + 3uv^)DI 
+ (4cru3 + Qavv^ + 12aviV2 + 8uu\ + Quvvi + 3uiw^) 

— 2uiD~^ o ui, 
A\^ = l2avDl + 20aviDl + {lQav2 + 12uv + v^)D^ 

+ 2(2crw3 + 2uvi + 3miw) — 2mi_D~^ o vi, 
A^^ = l2avDl + l&aviDl + {12av2 + 12uv + v^)D^ 

+ (4ctw3 + 8uvi + 6uiv + 3v^vi) — 2viD^^ o ui, 
A\^ = 8aDl + 2(4m + 3v^)D^ + 2{2ui + 3vvi) - 2viD-^ o vi. 

The matrix elements of A^ are given as 

^fi = 32a\Dl + 80a\iD^ + 8c7(14c7i;2 + 8uv + v^)DI 
+ 4a{22aV3 + 24uvi + 24wit; + 9v'^vi)Dl 

+ 20auv2 + 26auiVi + 16(TM2f + 9av^V2 + 12avvi 
+ 2uv^)Dx +4{2(7'^V5 + 6(TMf3 + llauiV2 + 9aU2Vi 
+ Aausv + 3av^V3 + 12(jvv\V2 + 3avl + Au^vi + 8uuiV 
+ 3uv'^vi + uiv^) 

— 4((TW3 + uvi + uiv)D~^ ou\ — 4:UiD~^ o {av^ + uvi + uiv), 
AI2 = IGa'^Dl + 8a(4u + 3v'^)Dl + 16a{3ui + 5vvi)Dl 

+ (32(7U2 + Mavv2 + Uavf + 16u^ + 2Auv^ + v^)D^ 
+ 4(2(T'U3 + 4(7VV3 + 6fTt;it;2 + 4uui + Auvvi + Suiv"^) 

— 4{(7V3 + uvi + uiv)D~^ ovi — 4uiD~^ o [ui + vvi), 
All = +16a'^Dl + 8(7(4u + 3v'^)DI + 16(t(3ui + Avvi)Dl 

+ (32(TU2 + 48(TDU2 + 28avl + Wu^ + 2Auv^ + v'^)D^ 

+ 4(2c7'U3 + 4(jvv3 + 8crviV2 + 4uui + 8uvvi + Suiv"^ + v^vi) 

— 4(^1 + vvi)D~^ oui— 4viD~^ o (av3 + uvi + uiv), 
AI2 = 32avDl + A8aviDl + 8{4av2 + 4.uv + v^)D^ 

+ 4(2(JW3 + Auvi + 4uiv + 3v'^vi) 

— 4(ui + vvi)D~^ ovi — 4:ViD~^ o (ui + vvi). 
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The matrix elements of are given as 

Al^ ^ -32(7^dI + 16cr2(-6u - 5v^)Dl + 80a'^{-3ui - 5vvi)D* 
+ 2(T(-160o-?i2 - 280cri;w2 - 204crw? - 481*^ - 80uv^ - 5v'^)Dl 
+ 4tT(— 60(7U3 — llOavv^ — 216(Tfii;2 — 72uui — 120miiwi 

— GOuiV^ - 15v^vi)Dl + 2{-48a^U4, - lOOa'^vVi - 244a^viV3 

— 168cr^uf — 96auu2 — 200cruvv2 — 136auvi — 68cruf 

— 260auivvi — S0au2V^ — 30av^V2 — dOcrv^Vi — 16u^ 

— AOu\^ - 5uv^)D^ 

+ (-16c7^U5 - 40cr^t;u5 - UOa'^viv^ - 208a'^V2V3 - 48auus 

— 120auvv3 — 232cruviV2 — 88auiU2 — 220auiVV2 — 156auiVi 

— 180au2VVi — iOauav^ — 20av^V3 — 120crv'^viV2 — QOavvl 

— 48u^ui — 80u^vvi — SOuuiv^ — 20uv^vi — 5uiv'^) 

+ 2(4cru3 + 6avv3 + 12crz;iW2 + 6uui + 6uvvi + 3uiV^)D~^ o m 
+ 8{(TV3 + uvi + uiv)D^^ o (av3 + uvi + uiv) 
+ 2uiD^^ o (4(T7i3 + 6avv3 + 12(7WiU2 + Quui + Quwi + 3uiv^), 
Al^ = -80cr^vDl - 22A(7^viDl + 40cr(-8cru2 - 4uv - v^)Dl 
+ 8cr(-30crw3 - 32?iwi - 30uiw - 2'i,v^vx)Dl 
+ (— 96(7^W4 — Vd2auv2 — 2h%au\V\ — \&}ou2V — Vo^av^v^ 

— 220avvl - SOu^v - AOuv^ - v^)D^ + 4(-4(7\5 - Uauvz 

— 22au\V2 — 18c7'U2i'i — IOctusI) — IQcrv^v^ — 3Qavv\V2 — Qcrv\ 

— 8u^vi — 2Quu\v — lOuv^vi — 5uiv^) 

+ 2(4cru3 + 6avv3 + 12av\V2 + 6uui + Quvvi + 3uiv^)D~^ o v\ 
+ 8{av3 + uv\ + uiv)D~^ o [ui + vvi) 
+ 2uiD~^ o (4c7i>3 + 6uvi + 6uiv + 3v'^vi), 
Al-^ = -80a\Dl - 176a\iD^ + 8a{-28av2 - 20uv - 5v^)DI 

+ 16cr(-llcrw3 - lAuvi - 15uiv - lQv^vi)Dl 

+ {-80a'^V4 - 160auv2 - 22AauiVi - 160au2V - 120av'^V2 

— 140(717?;^ - 80v^v - 4Quv^ - v^)D^ + {-I&ct'^v., - ■i8auv3 

— 88(TWi'i72 — 88o"U2Wi — 40(7w,3i; — AQav^v^ — l&{)(jvv\V2 — 3Qav\ 

— A:8u^vi - 80uuiv - 80tmi'^vi - 20?iiw^ - bv^Vi) 
+ 2(4(Tf3 + duvi + %uiv + 3v^vi)D^^ o ui 

+ 8(mi + vvi)D^^ o {av3 + uvi + uiv) 

+ 2viD~^ o (4(7113 + Qcrvv3 + 12av\V2 + &uui + 6uvvi + 3uiv'^), 
A\^ = -32a'^Dl + 16c7(-4u - ^v^)Dl + 48c7(-2«i - 5oti)£>^ 
+ 2(-32(TU2 - 8Qavv2 - 52avl - 16m^ - AQuv^ - f)V^)D^ 
+ 4(— 4c7U3 — IQCFVV3 — l&aviV2 — 8uui — 20uvvi 

— lOuiv^ - 5v^vi) 

+ 2(4cru3 + 6uvi + 6uiv + 3v'^vi)D~^ o vi 
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+ 8{ui + vvi)D^^ o (ui + wi) 

+ 2viD^^ o {Aav3 + 6uvi + 6uiv + Sv'^vi). 

Similar to the previous cases, we checked the conditions for Hamiltonianity and 
compatibihty of ah six Hamiltonian structures. It is also easy to check that all six 
structures are related by the recursion operator constructed for symmetries of the 
Boussinesq equation, see, e.g., [7]. 

7. The coupled KdV-mKdV system 

We shall now describe a Hamiltonian structure for the coupled KdV-mKdV 
system of the form 

Vt = -Vxxx + Sv^Vx + 3uVx + SUxV. (19) 

This system arises as the so-called bosonic limit of the = 2, a = 1 supersymmetric 
extension of the KdV equation ; integrability properties of this system (existence 
of a recursion operator) were studied in JB] . In [3] , by means of the prolongations 
structure techniques, a Lax pair for (jl9p was constructed. 

Denote the evolution equation corresponding to (IT^ by S'°° and choose for 
coordinates in S"^ the functions 

X, t, u = uo, V = vo,...,Uk, Vk,-.- 

Then the total derivative operators restricted to S""" are written in the form 

^ d ir-^f d d 

Dx = ^ + >[ Uk+l ^ h Vk+l 



k>0 



where 



/ = — U3 + 6uui — 3vv3 — 3viV2 + 3uiv^ + 6uvvi, 
g = — U3 + iv^vi + iuvi + iuiv 
are the functions at the right-hand side of ([T9|. 
The i*^- covering. The linearization operator restricted to S'°° is 

'^11 ^12' 



,^21 £22 

where 

4i =Dt + Dl~ (6u + iv'^)D,, - 6(mi -t- wi), 

tx2 = 3vD'^ + 3viD^ — {6uv — 3v2)Dx ~ &u\v — %uv\ -I- 3w3, 

h\ -ivDx - 3wi, 

I22 ^Dt + Dl^ 3(u + v^)Dx - (6twi + 3u). 
Consequently, the adjoint operator is 

■til ■'^21 

n* n* 
*12 *22 

where 



(21) 



(22) 
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£*2 = ~3vDl - 6viDl + 6{uv - V2)D^, 

e;2^~Dt^Dl + s{u + v^)D,. 

Following the general theory of Section [3l we now construct the £ J-covering for 
the equation S""" by introducing new odd variables p ^ po, q ^ qg, . . . ,pk, qu, ■ ■ ■ , 
Pk = D^{p), Qk — D^{q): that obey the equations 

Pi = -P3 + (6u + 3t)^)pi + 3wgi, (23) 

qt = Svp3 - 6vip2 + 6{uv - V2)pi - 173 + 3(u + v^)qi. (24) 

Solving the defining equations. We now introduce a vector-function of the form 

'F\ = fE^iFtP^ + Frq^)\ 
,GJ [X^(GfP^+G-q^J' 

where i^", F", G", are functions on S"^, and solve the equation 

'£11 ii2\ fF\ ^ ^25) 



£21 £22/ \G 

The operators £ij here are lifted to the £J-covering, which means that the total 
derivatives are now of the form 

where /' and g' are the right-hand sides of and (P^ . respectively. 
The following solution was obtained 

F = -p3 + 4mpi + 2uipo + 2vqi, 

G = 2i>pi + 2vipa + qi, 

to which there corresponds the operator 

(-DI+ AuD, + 2ui 2t;i?,\ 

2i,i^, + 2«i i?, j- ^^^^ 

T/ie Hamiltonianity test. We shall check now that the operator A presented by (j27p 
is Hamiltonian. The first property is obvious: evidently, A* = —A, i.e., yl is a skew- 
adjoint operator. 

To check the second property, we construct the bivector 

Wa = Fpo + Gqo = (-p3 + Aupi + 2uipo + 2wgi)po + (2wpi + 2vipo + qi)qo 

= P0P3 - 4upoPi - 2vpoqi + 2vpiqo + 2vipoqo - q^qi 

and verify condition ([5]), i.e., 

^r5WA5WA^SWA5WA\^^ ^28) 
\ 5u Sp 5v 5q J 



But 

SWa 
6u 

SWa 
Sv 

SWa 
Sp 

SWa 

~S^ 



-4pogi, 

2{p3 - Aupi - 2uiPq - 2wgi), 
2{-2vpi - 2vipo - qi) 
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and consequently 

5Wa 5Wa 5Wa 5Wa 
5u 6p Sv 5q 

i.e., (EHl) holds. 



~8popiP3 = Dx{-8pQpip2), 



= A£iX), (29) 



Existence of a Hamiltonian. Let us show that the KdV-mKdV system (fTO)) possesses 
a Hamiltonian, i.e., its right-hand side may be represented in the form 

where A is the Hamiltonian operator described above and X is the dx-component 
of a conservation law ij — X dx + T dt (the energy) . 

We computed directly several conservation laws of lower order and obtained 
the following results (for the sake of briefness, we omit the corresponding dt- 
components): 



m ■ 


X 




m ■■ 


X 




m ■ 


X 


= i(u^ + Uv'^ - VV2), 


V6 ■ 


X 


= I2u^ + 2Av?v'^ ~ &UU2 



6uv'^ — 30UVV2 — 3U2V^ — 8v^V2 + 6W4. 



Generating functions corresponding to these conservation laws, that is, vector- 
functions of the form 

are 

(fix = 0, Vi = 1; 

^2 ==1, -02 = 0; 

(fiQ = 6{6u^ + Suv^ — 2u2 + v'^ "06 = 12(4u^w -I- 2uv^ — 5uv2 — 5uiVi 

— '6VV2 — v\)] — 3U2U — 4:v'^V2 — 4:Vvl + V4). 

Applying A to £-{X), where X corresponds to 774, we see that holds. 



Theorem 2. The coupled KdV-mKdV system (1191) is Hamiltonian with respect to 
the the Hamiltonian operator (|27p and possesses the Hamiltonian X — ^{u^ + uv'^ — 
VV2)- The corresponding energy is given by the form 

rj — ijijJ^ + "^^"^ ^ VV2) dx + 2 ^^^"^ ^ 9u^v^ — 2uu2 + 3uv^ — lluw2 + uv\ 

+ u\ ~ uivvi — Au2V^ — 6v^V2 — Sv'^v'l + VV4 — + V2) dt. 

This structure is unique in the class of Hamiltonian structures independent of x 
and t and polynomial in Uk, Wfc. 

Proof of uniqueness. Let us first note that equation (|19p admits a scaling symmetry 
that allows to assign gradings to all variables x, t, Uk, and Vk'. 

\x\ = -1, \t\ = -3, |ufe| = /c + 2, \vk\^k+l; 

the grading of a monomial is the sum of gradings of the factors entering this mono- 
mial. In particular, |/| = 5, \g\ = 4. All constructions are in agreement with these 
gradings and we may restrict computations to homogeneous components. Since the 
grading of the expression B£{X) is \B\ + \£(X)\, we conclude that the grading of 
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the generating function E{X) is less than that of the right-hand side of p^ . This 
fact restricts the choice of possible Hamiltonians just to several ones and by a direct 
computation we find that the only possible solution is given by Theorem [2j □ 

Discussion: nonlocalities. In spite of the previous result, we have constructed an- 
other Hamiltonian operator for the system under consideration. This operator ex- 
ists in an appropriate nonlocal setting. First, we introduce a new nonlocal variable 
w defined by 



Wt = 3UV + — V2 

and corresponding to the conservation law rji (see Example [5]). 

In this nonlocal setting, it is possible to extend the £J,-covering by adding odd 
nonlocal variables ri, r2, defined by the relations 

ri,x = qovi +P0U1, 

ri,t = ~q2Vl +P2i-Ul ~ 3vVi) + qiV2 + Pi{u2 + 3ot2 - 3w^) 

+ qo{3uvi + 3uiv + 3w^wi — 113) 

+ Po{6uui + 6uvvi + 3uiv^ — U3 — Svv^ — 3viV2); 

r2,x = ^(Jocos(2ii;)'(;i - ^qoSm{2w)u 

- Po cos(w)(iui + vvi) + Pa sm{2w){uv - ^^2), 

r2,t = ^^92 cos(2u))ui + ^q2 sm{2w)u + ^p2 cos(2w)(ui - Wi) 
+ sm{2w){uv + V2) - qi cos(2w)(uw + ^1^2) 

- qi sin(2ii;)(iui -|- wi) 

115 

- pi cos{2w){uv^ + -U2 + -VV2 + -vl) 

5 1 1 

-I- Pi sin(2w)(-uui + -uiv - v'^Vi - -V3) 

+ —qa cos(2u')(5uui -t- uiv + v^vi — v^) 

3 1 11 

+ qo sin(2ii;)(--u^ - -uv'^ + -u2 + -VV2 + vj) 

3 1 

-I- Po cos(2w)(— 3u7ii — 6uvvi — i^^iv^ + -7/3 — v^vi 
3 5 

+ l^VVz + -V1V2) 

5 3 3 

-I- Po sin(2u))(3u^w -I- uv'^ — -:^uv2 — 3uiV\ — -^U2V — 2^^^2 

- 3vvl + iu4); 

rz,x = -^(7ocos(2w)u - igo sin(2w)t;i 

+ Po cos{2w){uv - ^V2) + Po sm{2w){^ui + vvi), 
rz,t = \q2 cos{2w)u + i(72 sin(2ti;)wi 
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+ ^P2 cos{2w){uv + V2) + ^P2 sm{2w){-ui + vvi) 

— qi cos(2w)(-wi + vvi) + qi sm{2w){uv — -^^2) 

5 1 1 

+ pi cos{2w){-uvi + -uiv - v^vi - -U3) 

115 

+ pi sin(2u))(uu^ + -U2 + -W2 + ^^i) 

3 1 11 

+ go cos(2w)(--u^ - -uv"^ + -U2 + -vv2 + vf) 

+ -qo sin(2w)(— 5uui — uiv — v'^Vi + v^) 

5 3 3 

+ po coa{2w){3u^v + uv^ — 2^^^ ~ "^^^^^ ~ 2^^^ ~ 2^^^^ 

- Zvv\ + iu4) 

3 1 3 5 

+ Po sin(2u))(3MMi + 'ouvvi + 2""!"^ ~ ^^^^ ~ 2^^^ ~ 2^^^^-'' 

In this extended setting, equation (j25p acquires a new solution of the form 
'F\ _ (Y.^{Ffp, + F-q,-rFfrif 

yGj [J:^iG^P^ + G:<l^+Grn), 

where 

F — (16uui + 12uvvi + 6uiv^ — 2uj, — Qvvy, — bviV2)po 
+ (16u^ + 12ui;^ - 8u2 - 17^2 - ^vl)pi 

— 3(4ui + bvvi)p2 - (8u + 5i;^)p3 +^5 

+ (2ui;i + buiv — W3)go + (Hw + 2t;^ — 4i;2)qi - 5wig2 — 4wg3 
+ 2r3(— 2 cos(2w)mw + cos(2w;)u2 — sin(2w)Mi — 2 sin(2w)wi;i) 
+ 2r2(cos(2w)Mi + 2 cos(2?ii)OTi — 2 sin(2w)MW + sin(2w)w2) 

- 3riui, 

G — {%uvi + Quiv + &v^vi — 2v^)pq + (lluw + 2v^ — &V2)pi 

- 7viP2 - 4:vp3 + {ui + 5vvi)qo + {2u + 5v^)qi - 53 

— 3wiri + 2(— cos(2w)i;i + sin(2w)M)r2 
+ 2(cos(2u;)u + sin(2w)i;i)r3. 

In the conventional matrix-operator form this solution looks as follows 

A' = L + N, 

where 

rfLu Li2\ AT ^ (Nil ^12 
\L21 L22)' \N2l N22 

are 2 x 2-matrix operators corresponding to the local and nonlocal parts of A' , 
respectively, and having the following entries: 

Lii =Dl- (8?i + bv^)Dl - 3(4ui + bvvi)Dl 

+ (16-u^ + 12™^ - 8-1*2 - 17i;u2 - Avl)D^^ 

+ 16uui + 12uvvi + 6uif^ — 2u3 — 6wt;3 — 5wif2, 
L12 = -"^vDl - bviDl + {lluv + 2v^ - Av2)Dx + 2uvi + buiv - V3, 
L21 = -ivDl - IviDl + (lluw + 2v^ - &V2)D,^ 
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+ 9uvi + 6uiv + Gv'^vi — 2v3, 



J22 



-DI + (2m + bv^)Dx + {ui + bvvi) 



and 



iV2i = -2,Yl,D-' o r^-o - 4^1^:1/?-! o Yl, - 4Yl,D-' o Y,% 



whereas 

are symmetries of the coupled KdV-mKdV system (see [6]) presented in the form 

Y"j = - cos(2w)(^Mi + vvi) + sm{2w){uv - ^^2), 
Yii = — cos(2i(;)?;i — - sin(2it;)w; 

Y^2 = cos{2w){uv — ^^2) + sm(2ii;)(iui + wi), 
Y^2 ^ ^2 cos(2w)u — - sm(2w)ui. 

Remark 14. Expressfons for the entries of the operator N were obtained as it was 
indicated in Remark [T^ 

As above, simulating the techniques developed for the local theory, we have 
checked that 

lA',A'j=0 and lA',Aj=0, 
i.e., A and A' are compatible. 

Remark 15. Though system ([T5| is Hamiltonian with respect to A' , there does not 
exist the corresponding Hamiltonian. This can be proved using the same techniques 
we used in the proof of Theorem[2l Nevertheless, the following facts are valid. Recall 
that (|19p possesses a recursion operator 6 . Denote this operator by R and note 
that our Hamiltonian operators are related to each other by means of this recursion 
operator, i.e., 

A' ^RoA. 

Moreover, in the same way one can construct a whole hierarchy of pairwise com- 
patible Hamiltonian structures. On the other hand, R generates a hierarchy of 
equations in which (1191) is the first term. Then A' is a Hamiltonian structure for 
all other equations of this hierarchy and these equations possess Hamiltonians with 
respect to A'. 

To make our exposition self-contained, we describe the recursion operator R in 
the Appendix. 
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Appendix: Recursion operator for the coupled KdV-mKdV system [B] 

To construct the recursion operator, it needs to extend the nonlocal setting 
introduced above. Namely, we add three new nonlocal variables wi, wn, and W12 
defined by 

Wi^t — + 3UV — U2 — 3VV2', 

wii^x — cos{2w)'Wiv + sin(2w)w^, 

wii,t — cos{2w){3wuv + wv^ — wv2 + uvi — uiv — v'^vi) 
+ sin(2w)(4Mw^ + — vv2 ~ w^); 

''^12,3: = cos(2w)u^ — sin(2u')ww, 

Wi2,t = cos(2w)(4uw^ — — 2vv2 + vf) 

+ am{2w){~3wuv — wv^ + wv2 — uvi + v'^vi) 
(see Example [S]). Then i? is a 2 x 2-matrix operator, 

i? — 1 -^^12 I 

\^2i ^22/ 

with the entries 

i?ii = Dl-Au-v^ 

- Y,^^,D-' o r^,2 + Yl2Dx' o - Iy-oD-\ 
R12 = 2vDl + viDx - 3uv + 2v2 

- Yl,D7/ o ri,2 + Y^,2D-^ o ri^^ + Yl^D-\ 
R21 = -2v 

Yl,D-' o ^^l, + Yl,D-' o - ^Yl,D-\ 
R22 = Dl-2u-v^ 

- Yl,D7/ o + Yl2D-^ o , + Y,\D-\ 

where Yi^O: ^i.ii ^1,2 are the same symmetries that enter the expression for the 
nonlocal Hamiltonian structure A' , 12,1 is another symmetry with the components 

Y^i = cos(2u')(— wiiui — 2i(;iiWi + 2wi2uv — W12V2) 
+ sm{2w){2wiiuv — W11V2 + 'W12U1 + 2wi2Wfi) 

— 2uvi — 3uiv — 2v'^vi + V3, 

Y21 — cos(2u')(it;iitii — Wi2u) — sin(2w)(wiiu + W12V1) 

- {ui + VVl), 

while 

i^r,i = -sin(2u;), 

V-i,! = 2(2sin(2i«)i;-wi2), 

V'1,2 = -cos(2w), 

i/'i.2 = 2(2 cos(2u;)i; + wu). 
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Note that tpi^i = (V'l dV'ii) and i/'i,2 = (^"1,2; "01,2) generating functions 

for (nonlocal) conservation laws. 
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